Theoretical investigation into the possibility of very large moments in Fei 6 N 2 
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We examine the mystery of the disputed high-magnetization a"-Fei6N2 phase, employing the 
Heyd-Scuseria-Ernzerhof screened hybrid functional method, perturbative many-body corrections 
through the GW approximation, and onsite Coulomb correlations through the GGA+U method. 
We present a first-principles computation of the effective on-site Coulomb interaction (Hubbard 
U) between localized 3d electrons employing the constrained random-phase approximation (cRPA), 
finding only somewhat stronger on-site correlations than in bcc Fe. We find that the hybrid func- 
tional method, the GW approximation, and the GGA+U method (using parameters computed from 
cRPA) yield an average spin moment of 2.9, 2.6 - 2.7, and 2.7 (ib per Fe, respectively. 



I. INTRODUCTION 

Though discovered in 1951 by Jackf^ cv"-Fei6N 2 (with 
crystal structure pictured in Figure [T]) first drew the at- 
tention of the magnetics community in 1972. It was 
then, 20 years later, that Kim and Takahashi 2 reported 
polycrystalline, mixed-phase Fe-N films with a saturation 
magnetization exceeding that of both cv-Fe and Co3sFe65 
(~ 2 x 10 6 A/m). However, it took another 20 years 
for the result to be reproduced (and, in fact, surpassed) 
by Sugita et aP3 Throughout the 1980s and '90s, other 
measurements of Fei6N2 thin films were reported that 
generally did not find this large magnetic moment 

Concurrently, density-functional theory (DFT) elec- 
tronic structure calculations were performedj^nl^l finding 
the moment per Fe ion to be modestly increased with re- 
spect to bulk bcc Fe but far short of the 3.5 hb reported 
by Sugita et al. It was showrP that LSDA+tfil cal- 
culations could yield an average moment comparable to 
that of some experiments (~ 2.8/ib per Fe), but the pa- 
rameters (U ~ 3.94, 1.0, and 1.34 eV on the Ad, 4e, and 
8h sites, respectively, with J = J7/10) were obtained via 
an embedded-cluster method with a small screening con- 
stant and were not calculated from first principles. Ad- 
ditionally, the J parameter is smaller than usually con- 
sidered appropriate for transition metals (typically one 
chooses either an atomic-like J of about 0.9 eV or else a 
more screened J of about 0.6-0.7 eV). 

Recently, further experimental evidence for the large 
magnetization has arisen,^ as well as a companion theo- 
retical paper^ reporting enlarged Fe moments achieved 
using LSDA+U (using U = 1.0 eV for the 4d site, 4.0 
eV for the 4e and 8h sites, and J = U/10). Ji et al. 
motivate their parameters by proposing that the Fe sites 
in the N-Fe octahedra form strongly correlated clusters 
in a metallic Fe environment, choosing a small U for the 
(within their model) more metallic 4d sites and a large U 
(chosen to be intermediate between that of FeO and Fe) 
for the 4e and 8h sites. They suggest that this model is 
supported by XMCD spectra that show additional fea- 




FIG. 1. (color online) Crystal structure of FeigN2, belonging 
to space group M/mmm. We use the PBE-relaxed structure 
in all calculations with a = 5.72 A, c = 6.29 A, x = 0.243, and 
z — 0.294. We will frequently refer to the three inequivalent 
Fe Wyckoff sites: 4d (red), 4e (green), and 8h (blue) (the N 
sites are black). 



tures at the Fe sites not seen in bcc Fe or other Fe-N 
phases 

In the present work, we perform an extensive search 
for the proposed large magnetization; we calculate the 
hyperfine field at the three Fe sites and compare with 
published Mossbauer spectra; we search for additional 
energy minima at moments away from the theoretical 
prediction as a function of tetragonal distortion; we 
apply the HSE06 hybrid-functional methocPSl and the 
GW approximation 2 ^ as implemented in VASF 22 to a - 
Fei6N2, testing the two methods on bcc Fe to ensure 
that any enhancement of the moment we obtain is gen- 
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uine. Further, we compute the effective on-site Coulomb 
interaction (Hubbard U) between localized 3d electrons 
employing the constrained random-phase approximation 
(cRPA^pRm (as implemented in the SPEJp3 extension of 
the FLEUR 2 '-I code), allowing us to provide for the first 
time first-principles predictions for the U and J param- 
eters. Finally, we present new PBE+lPH calculations 
using these parameters and discuss their implications for 
existing models. 



II. COMPUTATIONAL DETAILS 

The hyperfine field calculation and the study of the 
dependence of the total energy on cell moment (fixed 
spin moment or FSM) and tetragonal distortion was per- 
formed using the FPLO codeP^We implemented the full 
relativistic expression for the hyperfine field 
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as given in Ref. [29| and references therein into FPLC 
Note, that here we give the general pre-factor to acco- 
modate for proper units. \v) are the solutions of the 

Kohn-Sham-Dirac equation, while a = f ^ q ~\ , with 

a being the vector of the Pauli matrices. fi n is the direc- 
tion of the nuclear spin moment. The wave functions ty v 
are expanded in local atom centered orbitals in FPLD. The 
effective integrand leads to a damping factor for ma- 
trix elements between orbitals from different sites, which 
allows to introduce the approximation that only terms 
with orbitals belonging to the atom at which the nuclear 
spins sits will be taken into account. The scalar rela- 
tivistic hyperfine fields only contains the Fermi contact 
term, while the full relativistic version contains all terms 
(including Fermi contact, orbital and spin dipole- nuclear 
dipole), due to the intrinsic 4-componcnt formulation of 
Eq. (JlJ and the use of 4-spinors. A non-relativistic limit 
of this expression reveals all of the separate terms. The 
major contributions come form the s-orbitals, for whom 
only the Fermi contact term contributes. Although, the 
core states contribute a large amount to the hyperfine 
field it has been showrP 1 ' that valence contributions can 
be sizable. The core contribution depends on the spin 
polarization of the core wave functions including the ef- 
fects of the crystal exchange potential and hence should 
be influenced and scaled by the local spin moment. In 
FPLD the semi-core (Fe 3s, 3p) states are treated like va- 
lence states. For this reason we include the valence and 
semi-core contributions via the onsite approximation as 
explained above. 

The FSM calculation was carried out within the PBE 
approximation using a 8 x 8 x 8 /c-point mesh in a linear 
tetrahedron method with Blochl corrections. Our VASP 
PBE+U (using the fully- local double-counting term), 
HSE06 and GW calculations used a plane-wave cut-off 



of 400 eV (29.4 Ry or 5.42 a^ 1 ). The PBE+U (HSE06 
and GW) calculations used an 8 x 8 x 8 (6 x 6 x 6) T- 
centered Monkhorst-Pack k mesh (also using the tetra- 
hedron method with Blochl corrections), employing a 
smaller 3x3x3 mesh for the exact-exchange sums. All 
of our VASP calculations use the PAW 32 pseudopotentials 
of Kresse and JoubertJ^ and all VASP moments are cal- 
culated within a sphere of radius 1.3 A on the Fe sites. 

To calculate the Hubbard U parameter we employ 
the constrained random-phase approximation (cRPA)PS 
within the full-potential linearized augmented-planc- 
wave (FLAPW) method using maximally localized Wan- 
nier functions (MLWFs) ! 25 l 34 l The cRPA approach offers 
an efficient way to calculate the effective Coulomb in- 
teraction U and allows to determine individual Coulomb 
matrix elements, e.g., on-site, off-site, inter-orbital, intra- 
orbital, and exchange as well as their frequency depen- 
dence. We use the FLAPW method as implemented 
in the FLEUR code^ with the PBE exchange-correlation 
potential for ground-state calculations. A dense 16 x 
16 x 16 k-point grid is us ed. T he MLWFs are constructed 
with the Wannier90 cod d 35 1 36 1 . The effective Coulomb po- 
tential is calculated within the recently developed cRPA 
methocP^ implemented in the SPEX code^ (for further 
technical details see Ref s. l2~4"l [231 and [37]). We use a 
3x3x3 fc-point grid in the cRPA calculations. 

In all calculations we use the PBE-relaxed structure 
with a = 5.72 A and c = 6.29 A (except in the FSM sur- 
vey) and internal parameters x = 0.243 and z = 0.294. 
As a final note, we consider all employed electronic struc- 
ture schemes (VASP, FLEUR, FPLD) to be equivalent with 
respect to numerical accuracy at the level required in 
the present study. The use of three different packages is 
motivated by the different implementations available in 
these codes. 



III. RESULTS AND DISCUSSION 

A. Hyperfine Field 

The hyperfine field provides a picture of the local mag- 
netic structure that, unlike measurements of the satu- 
ration magnetization, does not require accurate estima- 
tion of the volume of a sample or its component phases. 
Mossbauer sp ectroscopy has been performed in many 
previous worka^^EHfil^ and the hyperfine field has been 
calculated 2112113 ^ from DFT. Our calculated B h r (found 
along with our calculated Fe moments in Table fij) agrees 
well with these past results; we find that the Fe sites with 
N nearest-neighbors exhibit approximately the same field 
(-23 and -22 T on the 4e and Sh sites), while B hf = -31 
T for the 4d sites. If we note, as previous authors haveP 
that DFT underestimates the hyperfine field by a sub- 
stantial, though nearly static, amount (~ 8 T in this 
case), then we also find reasonable agreement with some 
of the experimental reports. Particularly, we agree well 
with Refs. [51 [351 an d [351 Although our hyperfine fields 
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agree numerically with those of Moriya et a/.pSl they 
claimed that the largest hyperfine field was to be found 
in the 8h site, a claim that is difficult to reconcile with 
the predicted relative magnitudes of the moments and the 
similar environment of the 4e and 8h sites. We note, how- 
ever, that this assignment of the hyperfine fields agrees 
better with the moments in the recent LSDA+U study 
of Ji et al.\m We cannot offer any new explanation for 
Sugita et a/.'s larger 46 T fielcPnor the presence of only 
one Mossbauer sextet in their later single-phase sample. 4 



Site Spin Moment (/xb) 


B hf (T) (scal.-rel.) 


Bhf (T) (rcl.) 


4d 2.85 


-34 


-31 


4e 2.17 


-25 


-23 


8h 2.36 


-25 


-22 



TABLE I. Spin moments calculated within scalar-relativistic 
PBE and hyperfine fields (both fully- and scalar-relativistic) 
for each Fe site calculated within PBE using the FPLO code. 



B. Fixed Spin-Moment Survey 

Generally, expansion of the lattice may not be an ef- 
ficient means of increasing the magnetization of a mate- 
rial, as the enhancement of the spin moments may not 
outpace the increase in volume. However, it is known 
that fee Fe, while ordinarily nonmagnetic, enters a high 
spin state upon expansion of the cell volumeP^ Therefore, 
we have explored the energy landscape as a function of 
total (spin) cell moment and ^, allowing the former to 
range from 34 - 48 Hb (corresponding to average spin 
moments of 2.12 - 3.0 /is per Fe) and the latter from 1.0 
-1.5 (holding a fixed in one set of calculations and vol- 
ume fixed in another). We only constrain the total spin 
moment of the cell and not the magnitude of the indi- 
vidual moments. In principle, the moments of the three 
inequivalent Fe sites could be arranged in many ways to 
obtain the same total spin moment; however, we simply 
accept the converged result for each structure and total 
moment without seeking out other possible minima. 

The results may be seen in Figure [2j We note that 
no additional local energy minima were observed apart 
from the PBE-relaxed structure (a = 5.72 A, c = 6.29 A, 
c/a Ri 1.10) and moment (2.44 /is). Although the energy 
minimum does tend to shift to higher moments with the 
increase of the volume through - , the enhancement is not 
sufficient to produce an increase in the magnetization. 
With a held fixed at a oxpt = 5.72 A, the average spin 
moment per Fe reaches 2.81 [is at - — 1.5, giving a 
magnetization of 1.49 x 10 6 A/m, compared to 1.77 x 10 6 
A/m at the experimental - and 1.75 x 10 6 A/m in bec Fe. 
If the volume is held fixed, the average moment does not 
depend strongly on ~ , remaining close to the PBE value 
throughout and decreasing to about 2.25fiB at ~ = L5. 
This supports the standard understanding of the LSDA- 



or GGA-predicted increase in the moment as arising from 
increased cell volume. 
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FIG. 2. (color online) Energy landscape of FS16N2 as a func- 
tion of tetragonal distortion and average moment per Fe with 
a) a held fixed and b) volume held fixed. Each contour rep- 
resents an increase of a) 10 mHa (0.27 eV) or b) 5 mHa 
(0.14 eV). There are no additional local energy minima in 
the parameter space examined. For fixed a, - = 1.5 gives 
l-iavg = 2.81/is- Although the average moment is higher at 
this point, it does not overcome the increase in volume, and 
the magnetization is only 84% of the magnetization of the 
PBE structure. For fixed volume, the average moment re- 
mains close to the PBE value, decreasing slightly as c in- 
creases above ~ 1.1a. 



C. HSE06 and GW 

It is possible that DFT cannot fully account for the 
physics that would give rise to greatly enhanced magne- 
tization in a"-Fei6N 2 , so we have also considered meth- 
ods that have arisen since the last wave of theoretical 
investigation into this material subsided. The HSE06 
screened hybrid functional method entails only a moder- 
ate increase in computational time with respect to PBE, 
and the inclusion of a static screening parameter for the 
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exact exchange term allows for the treatment of metal- 
lic systems — unlike the parent Hartree-Fock method — as 
well as speeding up calculation further. HSE06 follows 
PBEO in its formulation of the exchange-correlation en- 
ergy, given by 
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The aforementioned screening parameter /i — 0.2 A -1 
partitions the exchange term into a short-range and a 
long-range component, achieved by appending erfc(//r) 
(the complementary error function) to the short-range 
terms and erf(/xr) to the long-range termP^ 

The GW approximation improves upon Hartree-Fock 
by treating electrons as dressed quasiparticles interacting 
via a screened Coulomb operator W . This replaces the 
purely real exchange-correlation potential with a complex 
self energy E = —iGW. In the initial step, the Green's 
function G and the screened Coulomb operator W are 
calculated from the wave functions obtained from a con- 
verged DFT calculation. The computation of W via the 
RPA is time-consuming, and consequently sonic short- 
cuts are sometimes employed. So-called "one-shot" GW 
or Go Wo is performed by calculating the quasiparticle 
energies using only these initial q uantit ies and yields im- 
proved results compared to LSDAP^EIl Nevertheless, the 
"one-shot" method still underestimates band-gaps due to 
the inaccuracies inherent in using an LSDA-obtained W, 
and improvement can be obtained by iterating G and W 
to self-consistency. We present results from G W , GW , 
and GW in this work. 

Figure [3] shows the partial density of states (pDOS) 
of each Fe site in HSE06 and GW. For comparison, 
we include the PBE-calculated pDOS for FeigN 2 and a 
fictitious "FeigNo" structure obtained by removing the 
N atoms without relaxing the structure. This latter 
case shows that, within PBE, Fe approaches the strong 
ferromagnetic state, with the majority d states nearly 
fully occupied, upon the N-induced volume expansion, 
yielding an average moment of 2.56 (Ib per Fe and a 
magnetization of 1.84 x 10 6 A/m, about a 5% increase 
over bcc Fe (with a bulk magnetization of 1.75 x 10 6 
A/m). The HSE06 pDOS shows a greatly enhanced ex- 
change splitting with respect to PBE-Fei6No and GW- 
Fei 6 N 2 , leading to an average moment of 2.86 fig per Fe 
(M = 2.06 x 10 6 A/m), whereas GW yields a more mod- 
erate 2.57-2.70 fi B per Fe (M = 1.85-1.95 x 10 6 A/m). 
The calculated spin moment at each site can be found in 
Table HI 

In the absence of experimental photo- or x-ray- 
emission data to which to compare, we must test the 
validity of the calculated moments by calculating the 
moments of better established materials. The last col- 
umn in this table shows the calculated spin moment for 
bcc Fe from PBE, PBE+U (which will be discussed in 
detail in the following section), HSE06, and GW. Our 
HSE06 result for bcc Fe agrees with previous worlsPD anc l 
demonstrates that, although the screened hybrid func- 



tional method improves on the Hartree-Fock treatment of 
metallic systems, it can overestimate the strength of the 
exchange and yield un-physical high spin states. How- 
ever, we also note that the calculated bcc Fe spin moment 
is not necessarily directly proportional to the calculated 
moments in a"-Fei6N2, so it is possible that the bcc Fe 
moment does not completely determine the accuracy of 
a method in this case. 



D. cRPA and PBE+U 

Previous attempts to explain the experiments that find 
high-magnetization have turned to LDA+U to describe 
the correlation effects that may be present in a"-FeigN2. 
However, as no first-principles calculations of the interac- 
tion parameters existed, it was necessary to motivate the 
choice of U (and J) by analogy with other systems or by 
applying a model. In particular, the explanati on for the 
enhanced magnetization proposed by Ji et and 
Wang et aZ™ requires that the Fe sites with N nearest- 
neighbors be more strongly-correlated than the Ad sites, 
which have no N neighbors. Without a set of firmly- 
established parameters, it is difficult to progress in un- 
derstanding this system, as the calculated moment is di- 
rectly dependent on U and J (see, e.g. Figure 3 in Ref. 
□U). 

Recently, the cRPA has been proposed as a first- 
principles method of obtaining the screened Coulomb 
matrix within a Wannier basis E2H2S1 Within the RPA, 
the polarizability P can be written 
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where the ipi an d £i are the DFT wave functions and 
their eigenvalues, and a runs over both spin channels. 
If one separates P into Pd, containing the correlated or- 
bitals, and P r , containing the rest, and if one co nsider s 
the unscreened Coulomb operator v, one can writ c! 23 * 2 ^ 



U = [1 - vP r ]- x v 

u = [i- up d ]- x u 



(4) 
(5) 



The matrix elements of the effective Coulomb potential 
U in the MLWF basis are given by 

U Rnin3]nin2 {ui) = J J < iR (r)w„ 3R (r)[/(r,r';w) 

x< 4R (r>n 2 R(rV 3 ^V, (6) 

where w n ji(r) is the MLWF at site R with orbital in- 
dex n and t/(r,r';w) is calculated within the cRPA. 
Strictly speaking, the Wannier functions are spin depen- 
dent. However, we find that this spin dependence affects 
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Method Ad Site (/is) 4e Site (/is) 8ft Site (/is) Average (/is) bcc Fe (/is) 

PBE 2^84 2T9 2~38 2^44 2^23 

PBE+U 3.08 2.62 2.74 2.71 2.67 

HSE06 3.06 2.83 2.91 2.86 2.85 

G Wo 2.90 2.31 2.49 2.57 2.33 

GW 2.95 2.35 2.53 2.64 2.62 

GW 2.96 2.41 2.57 2.66 2.65 

GW(s,pval.) 3.00 2.50 2.64 2.70 2.59 



TABLE II. Calculated spin moments for all methods presented in this work. The PBE+U results were obtained using the 
cRPA-obtained interaction parameters U = 3.99, 3.12, and 3.52 eV for the 4d, 4e, and 8ft sites, respectively (J = 0.64, 0.59, 
and 0.61 eV), except for bcc Fe, for which we used U — 3.16 eV and J = 0.68 eV as in Ref. 1251 GW (s,p val) denotes a 
VASP-GW calculation in which the Fe 3s and 3p electrons are treated on the same level as the 3d and 4s. 



the values only little. For simplicity, we ignore the spin 
dependence here and give the spin-averaged values in the 
following. 

In our Spex-cRPA calculation, we choose the Fe d or- 
bitals as our correlated subspace and compute the in- 
teraction parameters found in Table [Tll*l Quantities with 
tildes are obtained from the fully screened Coulomb ma- 
trix U, while plain symbols are the sp-screened quantities 
that enter into the PBE+U calculations. The U, U' , and 
J (and their fully-screened counterparts) are averaged at 
each site as follows: 



Udft+u = F° = — ^2 U, 

m 
m < n 
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(7b) 
(7c) 
(7d) 



Site 


Udft+u 


U 


U' 


J 


U 


U' 


J 


4d 


3.99 


5.02 


3.74 


0.64 


1.80 


0.71 


0.53 


4e 


3.12 


4.14 


2.95 


0.59 


1.56 


0.55 


0.49 


8ft 


3.52 


4.50 


3.27 


0.61 


1.68 


0.62 


0.51 



TABLE III. The calculated on-site interaction parameters (all 
in eV) from cRPA for a"-Fei6N2, showing a small increase in 
correlation with respect to bcc Fe. Quantities with a tilde are 
computed from the fully-screened Coulomb potential, while 
plain quantities are computed from the partially screened po- 
tential (omitting d screening). Udft+u is the U parameter 
that enters into the PBE+U calculations. 



Ref. [25]— U = 3.16 eV and J = 0.68 eV. We use the fully- 
local (FLL) double counting correction in the calculation 
of both the bcc Fe and the Fei6N2 moments. Although 
this choice may seem strange in metallic systems, the 
around-mean-field (AMF) term opposes the formation of 
moments in generate and here produces moments ~ l/i^ 
below the expected value in bcc Fe. It should be noted, 
however, that the choice of the double counting term in 
PBE + U is not unique and thus leaves an ambiguity in 
the calculated moments even if U and J were computed 
with a well-defined method.. 



E. Orbital Moment 

In solids, the orbital moment is typically nearly 
quenched, but in some extreme cases, such as UN, 47 
the orbital moment can be comparable to the spin mo- 
ment. PBE calculations give an orbital moment per Fe 
of only 0.05 hb in bcc Fe (Table IV I, but this may be 
increased somewhat in Fei6N2. To explore this possi- 
bility, we calculated the orbital moment within PBE, 
PBE+orbital polarization correction (OPC)J^ PBE+U 
(using the cRPA parameters), and "one-shot" Go Wo us- 
ing FPL0 (for the OPC calculation) and VASP (for the 
rest). Each method shows a small increase in orbital 
moment compared to bcc Fe, yielding about 0.1 - 0.2 
fiB per Fe atom and an increase of 0.01 - 0.05 /ib over 
bcc Fe. This small increase cannot explain those results 
that claim average Fe moments in excess of 3 /i#. Our 
PBE+U and Go Wo results predict average total (spin + 
orbital) moments of 2.88 and 2.63 [Ib, respectively. 



IV. SUMMARY 



We note that these parameters differ both quantita- 
tively and qualitatively from previously proposed mod- 
els, particularly those that suggest large differences in 
correlation strength between Fe sites. The spin moments 
from PBE+U, for FeigN2 as well as bcc Fe, can be found 
in Table [TT] The PBE+U spin moment for bcc Fe was 
calculated using the interaction parameters computed in 



We have examined the electronic and magnetic struc- 
ture of a"-Fe 16 N 2 within PBE, PBE+U, HSE06, and 
GW. Within PBE, we find spin moments and hypcrfine 
fields that agree with past results, and we do not find that 
any high-magnetization state arises as ~ changes from the 
experimental value. We have provided effective Coulomb 
interaction parameters calculated via cRPA and have 



6 



Method 


4*7 Sitp ( a id) 


Ap fti+p ( 1 1 


Rh Site in™} 


A ^/"pra cfp ( 1 1 t>) 
nvcidgc \H J £> ) 


Hpf 1 ill T~> 1 


PBE 


0.06 


0.06 


0.05 


0.06 


0.05 


PBE+U 


0.20 


0.16 


0.16 


0.17 


0.12 


PBE+OPC 


0.09 


0.11 


0.10 


0.10 


0.09 


GoWo 


0.06 


0.06 


0.05 


0.06 


0.05 



TABLE IV. Calculated orbital moments in Fei 6 N 2 within PBE, PBE+U, PBE+OPC, and G W . The orbital moment is 
increased by only 0.01 - 0.03 [ib per Fe with respect to bcc Fe. 



used them in our PBE+U calculations. We find that 
PBE+U and HSE06 gives average spin moments per Fe 
of 2.71 and 2.86 [Ib but also greatly overestimate the mo- 
ment of bcc Fe (experimentally about 2.2 ^b)- GW gives 
smaller moments, 2.57 - 2.70 [is per Fe, a slight increase 
over the PBE moment. G W , GW , and GW all over- 
estimate the bcc Fe spin moment by different amounts 
despite their similar predictions for Fei6N2, with Go Wo 
giving the most reasonable bcc Fe moment due to its 
close dependence on the PBE result. In all cases, we find 
that the 4e and 8h sites have smaller moments than that 
on the Ad site. 

We have also presented calculations of the orbital mo- 
ment on the Fe sites obtained within PBE, PBE+OPC, 
PBE+U, and Go Wo- We find that the orbital moment 
is not completely quenched and may add 0.1 - 0.2 /is to 
the average total moment per Fe, a small increase over 
bcc Fe. 

In order to evaluate the varying results found above, 
one must understand the purposes of and approxima- 
tions inherent in the methods presented. In addition to 
the shortcomings of the mean-field-like treatment of cor- 
relations within PBE+U, there are two notable avenues 
for error in this method: the need to choose the U and J 
parameters and the lack of a priori justification for the 
double-counting corrections. Dependence on the choice 
of interaction parameters is not a fundamental problem 
and can be alleviated as we have done here by computing 
them through some appropriate first-principles method. 
The choice between the FLL or AMF double-counting 
corrections, while straightforward when treating insula- 



tors, can be less obvious in semi-localized magnetic sys- 
tems, and furthermore no method exists for determining 
the exact form of the correction. The hybrid functional 
method's dependence on parameters is fundamental to 
the approach, although it is mitigated somewhat by the 
use of predetermined parameters such as in HSE06. How- 
ever, these parameters were primarily chosen to produce 
reasonable band gaps and may need to be altered to prop- 
erly treat metallic systems. In principle, the GW approx- 
imation should be the most accurate of those presented 
here. The Go Wo and GWo methods maintain good con- 
tact with the PBE results while incorporating first-order 
exchange and correlation effects. However, some care 
must still be taken; we have shown that the results do 
depend on which electrons are treated as valence and 
which are absorbed into the core pseudopotcntial. Lastly, 
we note the need for additional, repeatable experiments 
that probe the electronic structure of the material in or- 
der to provide a better basis for comparison with theory. 
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